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ABSTRACT 



Context. Among multi-planet planetary systems there are a large fraction of resonant systems. Studying the dynamics and formation 
of these systems can provide valuable informations on processes taking place in protoplanetary disks where the planets are thought 
have been formed. The recently discovered resonant system HD 60532 is the only confirmed case, in which the central star hosts a 
pair of giant planets in 3:1 mean motion resonance. 

Aims. We intend to provide a physical scenario for the formation of HD 60532, which is consistent with the orbital solutions derived 
from the radial velocity measurements. Observations indicate that the system is in an antisymmetric configuration, while previous 
theroretical investigations indicate an asymmetric equilibrium state. The paper aims at answering this discrepancy as well. 
Methods. We performed two-dimensional hydrodynamical simulations of thin disks with an embedded pair of massive planets. 
Additionally, migration and resonant capture are studied by gravitational N-body simulations that apply properly parametrized non- 
conservative forces. 

Results. Our simulations suggest that the capture into the 3:1 mean motion resonance takes place only for higher planetary masses, 
thus favouring orbital solutions having relatively smaller inclination (z = 20°). The system formed by numerical simulations qual- 
itatively show the same behaviour as HD 60532. We also find that the presence of an inner disk (between the inner planet and the 
star) plays a very important role in determining the final configurations of resonant planetary systems. Its damping effect on the inner 
planet's eccentricity is responsible for the observed antisymmetric state of HD 60532. 

Key words, planets and satellites: formation, celestial mechanics, hydrodynamics, methods: N-body simulations 



1. Introduction 

A significant fraction of multi-planet planetary systems contain 
a pair of giant planets engaged in a mean motion resonance 
(MMR). These planets are mainly in the low order 2:1 MMR 
(as GL 876, HD 128311, and HD 73526), but in a few systems 
higher order resonances have been suggested as well; two plan- 
ets around 55 Cancri may be in 3:1 MMR, or HD 202206 may 
host a pair of planets in 5:1 MMR. In the recently discovered 
system, HD 45364, the giant planets are revolving in 3:2 MMR 
ICorreia et alJ (120091) . 

The observed orbital solutions and the formation of the ma- 
jority of resonant systems have been studied t horoughly by many 
authors. It has been shown, for instance, by iKlev et all (12004]) 
that a sufficiently slow migration process of two giant planets 
embedded in an ambient protoplanetary accretion disk ends with 
either a 3: 1 or 2: 1 resonant configuration depending on the speed 
of migration. The formation of resonant systems in 2:1 MMR 
has been modelled exhaustively by hydrodynamical and N-body 
simul ations as well. T he system GJ 876 h as been inve s tigate d 
by iLee & Peald d2002l) . Kley et al (2005), ICrida et all 12 008), 
and th e sy stems HD 128311 a nd HD 73526 by San dor & Klevl 
(2006) and lSandor et al.l (120071) . respectively. Most recently, the 
formation of the system HP 45364 w i th plan ets in the 3:2 MMR 
has been investigated by iRein et al.l (|2010). Analytical studies 
related to the stationary solutions of th e 2:1 and 3:1 M MR have 
been done bv lBeauge etal] d2003l) and lBeauge et alJ d2006l) . for 
instance. 



Regarding the 3:1 MMR case, the pres ence of the plane t 
55 Cancri-c has already been questioned dNaef et alJ [2004). 
Additionally, recent orbital fits indicate that the planets 
55 C ancri-c and 55 Ca ncri-d are not in a resonant configura- 
tion dFischer et alJ l2008). Thus, until the recent discovery of the 
two planets in a 3: 1 MMR around the F-type star HD 60532 by 
Deso rt et al ] d2008h . there has been a serious lack of knowledge 
about the observed behaviour of 3:1 reson ant systems. The re- 
sults of the dynamical study performed by [Desort et alJ {2008) 
did not prove without any doubt that the giant planets are in fact 
in a 3:1 MMR. 

The final confirmation of the 3:1 MMR between the gi- 
ant planets in HD 60 532 is given in a recent paper by 
lLaskar & Correial (|2009), in which two new orbita l fits are pro- 
vided, slightly improving on the previous fit of iDesort et all 
(2008) with i = 90°. Th rough a detailed s t ability analysis of 
different orbit integrations. lLaskar & Correial d2009l) suggest that 
an inclination of i = 20° is the most likely configuration for a co- 
planar model. By this assumption, the planetary m asses are in- 
crease d by a factor of 1 / si n(Q a 3 in compar i son tolDesort et all 
(2008}). In a recent paper. IHbert & Tsiganisl d2009l) analyse the 
excitation of mutual inclination for planetary systems driven into 
resonance by planet-disk interaction. They use a damped N- 
body evolution and find that for stronger eccentricity damping, 
the mutual inclination is less excited. However, they could not 
place any constraints on the observed inclination of the system 
HD 60532. 

The very large difference in planetary masses between the 
two orbital solutions having i = 20° or 90° also raises the impor- 
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Table 1. Orbital solutions provided bv lLaskar & Correial J2009) 
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tant question of which orbital solutions are preferred by forma- 
tion based on the planetary migration scenario. In the present pa- 
per we aim at answering this question by performing fully hydro- 
dynamical evolution of planets embedded in the disk. Through 
this procedure, we obtain realistic migration and eccentricity 
damping rates that will allow us to determine the most proba- 
ble final state of the system. 

The paper is organized as follows. First, we numerically in- 
tegrate the orbits o f giant planets u s ing th e two sets of orbital 
solutions given bv lLaskar & Correial d2009l) as initial conditions. 
Then by considering different planetary masses, we investigate 
the possible capture into the 3:1 MMR between the giant plan- 
ets. After having formed a resonant system, we compare the re- 
sults of our simulations to the orbital behaviour of giant planets 
obtained from numerical integrations of Sect. 2. Finally, by per- 
forming gravitational three-body numerical integration with dis- 
sipative forces for migration, we study how an inner disk influ- 
ences the behaviour of the system toward its stationary solutions. 
We demonstrate that the presence of the inner disk determines 
the final resonant configuration of the system. 



2. Orbital behaviour of the system HD 60532 

In this section we present the orbital behaviour of the giant plan- 
ets around the star HD 60532. To integrate the equations of the 
three-body problem, we used a Bulirsch-Stoer numerical inte- 
grator. The initial conditions for o ur simulations were calcu lated 
from the orbital solutions given bv lLaskar & Correial (120091) . see 
Fits I and II in Table Q] 

A resonant system can be characterized by studying the be- 
haviou r of the resonant angles. According to iMurrav & Dermottl 
( 1999), in the planar case of the 3 : 1 MMR, there are three dif- 
ferent resonant angles: 8\ - r iA2-A\-co\ -ci>2, #2 = 3A-2-A\-2u>2, 
and 63 = 3A2 - Ai — 2wi, where /I, denotes the mean longitude 
of the z'-th planet, and to, its periastron. Here index "1" stands for 
the inner and index "2" for the outer planet. By introducing the 
corotational angle Aco - a>2 - co\ and arbitrarily choosing a 8,, 
all the remaining resonant angles can be expressed with 0, and 
Au>. If one of 0, (z = 1 , 2, 3) librates around a constant value, the 
system is in a 3 : 1 MMR. If Aco also librates, the system is said 
to be in apsidal corotation. 

By studying Fig. [TJ one can follow the behaviour of the 
eccentricities and resonant angles as functions of time for Fit 
II. The eccentricities show quite large oscillations, the inner 
planet's eccentricity varies between e\ ~ 0.1 - 0.33, while the 
outer planet's eccentricity is between e2 ~ 0-0.15. All reso- 
nant angles are librating, 8\ around 0°, 62, 83, and Alo around 
180°. This means that the system is in 3:1 MMR and also in ap- 
sidal corotation, and furthermore the orbits of the giant planets 
are anti-aligned. In Fig. [1] we display the behavio ur of the reso- 
nant angle 83, (as used bv lLaskar & Co rreia 2009) and Alo. By 
using Fit I , we found the same dyn amical behaviour of the giant 
planets as lLaskar & Correial d2009l) . therefore we do not display 
these results here. 
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Fig. 1. Dynamical behaviour of the system HD 60532, where 
the initial con ditions of the numerical integration are calculated 
from Fit II of lLaskar & Correial 42009). The upper panel shows 
the evolution of the eccentricities (inner planet red curve, outer 
planet green curve), middle panel the behaviour of the resonant 
angle #3, while the bottom panel shows the time evolution of the 
corotation angle Alo. 



Regarding the behaviour of the resonant angles and the pos- 
sible apsidal corotation, a co mprehensive analytic al study of the 
3:1 MMR has been done by iBeauge et al.1 (120031) . The authors 
find that the resonant angle 83 and Aco both librate around 180° 
if the eccentricity of the inner planet is below a certain limit: 
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e\ ss 0.13. Above this limit (i.e. e\ > 0.13), there appears asym- 
metric libration of the resonant angles, whose mean values de- 
pend on the eccentr icity of the inner planet. This result has also 
been confirmed by Kiev et aD d2004l) . where the mean value of 
the corotation angle was found to be Aa> ss 110°. Returning to 
the behaviour of the system HD 60532 (obtained from radial 
velocity measurements), although the inner planet's eccentric- 
ity is higher than the limit e\ = 0.13, the corotation angle does 
not yield any asymmetric libration; on the contrary, it librates 
around 180°, see the bottom panel of Fig. Q] The reason for this 
behaviour might be that, for some reas on, the system coul d not 
reach the stationary solution found by Beauge et al. (2003) dur- 
ing the planetary migration. A possible explanation of this phe- 
nomenon will be presented in Sect. 5. 

In what follows, we attempt to reproduce the dynamical be- 
haviour of the system shown in Fig.[T]by a migratory formation 
scenario using hydrodynamical and gravitational 3-body simu- 
lations with properly parametrized non-conservative forces for 
migration. 

3. Numerical setup of hydrodynamical simulations 

To study the formation of the resonant system HD 60532 through 
migration of the giant planets, we performed a series of full 
hydrodynamical simulations of a protoplanetary accretion disk 
containing two embedded giant planets. For the hydrodynamical 
simulations we used the locall y isothermal v ersion of the FARGO 
public hydrodynamical code (Masset 2000) and the code RH2D 
dKlevlll999l) both well-suited to investigating planet-disk inter- 
actions. 

For our simulations, we used a flat accretion disk with an 
aspect ratio H/r = 0.05 extending between r m ,„ = 0.2 and 
fmax - 5.0 in dimensionless units. The disk has an initial surface 
density profile S(r) = l,or ~ 1 ^ 2 , where the dimensionless surface 
density value at one distance unit is £q = 3.1 x 10~ 4 . This value 
corresponds to 1% of the stellar mass (M» = 1.44M ). During 
our simulations, the accretion of the gas in the protoplanetary 
disk was driven by an alpha-type viscosity with a = 0.01. A 
torque cut-off around the planet of Q.6RR OC i ie was also applied. 
The computational domain (bounded by r„„„ and r max ) was cov- 
ered by 256 radial and 500 azimuthal gridcells. The radial spac- 
ing is logarithmic, while the azimuthal is equidistant, resulting in 
nearly quadratic gridcells (meaning that at each radius rA(f> ~ Ar, 
approximately). 

The planets were initially placed at r = 1.0 and 2.5 dis- 
tance units from the star. To avoid the transient effects and ob- 
tain a quasi steady state initial setup, we first integrated the disk- 
planets-star system keeping the orbits of giant planets fixed on 
circular orbits until 500 orbital periods of the inner planet. The 
surface density distribution of the accretion disk obtained after 
the above integration time is shown in Fig. [2] where the plan- 
etary masses correspond to the higher masses of Fit II, and we 
used a logarithmic scale for the surface density values for better 
visibility. 

There is an essential differ ence between our hydrodynami- 
cal setup and the one used by iKlev et al.l (2004) (Model C in 
the cited paper). In the latter model, the initial surface density 
profile tends to zero with r — » 0, because the inner disk was as- 
sumed to be cleared by the accretion of the gas to the star. In 
such a case, the outer planet feels the negative torques of only 
the outer disk and migrates inward until it captures the inner 
planet into the 3:1 MMR, which revolves in a practically gas- 
fre e environment. On t he o ther hand, a s demonstrated recently 
bv lSandoretai1(l2007h and lCrida etaP (120081) . the presence of 
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Fig. 2. The surface density distribution of the protoplanetary ac- 
cretion disk after 500 periods of the inner giant planet. The disk 
model is described in Sect. 3, the planetary masses correspond to 
Fit II. Light colours are for higher, while dark colour shades are 
for low values of the surface density. For a better visualization 
we used a logarithmic scale. 

the inner disk plays a very important role in damping the eccen- 
tricity of the inner planet during the migration process, hence 
in shaping the final behaviour of the system. Therefore in our 
disk model we did not allow the inner disk to be emptied, so the 
planets (mainly the inner one) can feel its effect. 

Having obtained almost cleared gaps after 500 orbital revo- 
lutions of the inner planet, the planets were released and their 
orbital evolution was followed until 4000 periods of the inner 
pl anet. For the planetary m asses we used the values published 
bv lLaskar & Correia d2009l) (see Fits I and II), and also an inter- 
mediate value calculated for i = 30°. 

4. Resonant capture for different masses of the 
giant planets 

In this section we show the results of our hydrodynamical simu- 
lations obtained by using different planetary masses. 

First we investigated the behaviour of the giant planets hav- 
ing low masses, corresponding to the Fit I, where i = 90°. 
In this case m\ = 1.048My, ni2 = 2.487My corresponding to 
qx = 0.0007 and q2 = 0.00165 dimensionless mass units (the 
stellar mass is unity). We used index "1" for the inner and index 
"2" for the outer planet. After releasing them, however, the plan- 
ets did not show a convergent migration, and no resonant capture 
took place. This was also true for the intermediate value of the 
inclination i = 30°, where the planetary masses are a factor of 
two higher than in Fit I. 

Finally, we used planetary masses corresponding to Fit II, 
where i = 20°, and the planetary masses are m.\ = 3.15My and 
in 2 = 7.46My or, in dimensionless mass units, q\ = 0.0021 and 
q2 - 0.005. (These values are quite close to the planetary masses 
used by Kley et al. 2004.) In this case the giant planets also be- 
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Fig. 3. The azimuthally averaged surface density profile of the 
protoplanetary disk after 500 periods of the inner giant planet for 
different planetary masses. Dashed line corresponds to planetary 
masses of Fit I, dotted line for the intermediate values, which are 
twice the masses of Fit I, and solid line for masses of Fit II. 



gin to migrate, but contrary to the previous two cases, the outer 
planet migrates faster than the inner one resulting in a resonant 
capture. 

Before describing the dynamical behaviour of the system af- 
ter the capture into the resonance, we should comment on why 
this phenomenon does not occur in the cases of lower mass plan- 
ets. We recall that for a resonant capture the planets should ex- 
hibit convergent migration. Consequently, if the planets migrate 
inward with the same speed, or if the inner planet migrates even 
faster than the outer one, no resonant capture will occur. For 
example, if the inner planet does not open a sufficiently empty 
gap, it may migrate faster than a usual type II migration. It is 
also true that the inner planet's inward migration is mainly gov- 
erned by the torques coming from the middle part of the disk, 
which is between the planets. (The outer disk does not really 
exert torques to the inner planet because it is separated by the 
wide gap opened by the outer planet.) On the other hand, the 
inward migration of the inner planet can be s topped and even re- 
versed by a sufficiently massive inner disk ( San dor etalJl2007t 
ICrida et al.ll2008h . From our hydro simulations one can see that 
only the massive inner planet of Fit II can open a sufficiently 
deep gap for a slower (type II) migration, and could also cre- 
ate a massive inner disk as pushing the disk's material toward 
the star. The planets with masses from Fit I and the intermedi- 
ate case (with i - 30°) open a shallower gap, the middle disk 
is more massive, and the inner disk is less massive for planets 
having masses from Fit II. This can be checked in Fig. [3] where 
we display the azimuthally averaged surface density profiles of 
the disk for the different planetary masses. Thus, as can be seen 
in the upper panel of Fig. [4] the migration of the inner planet 
(of Fit II) is slow in the beginning and is stopped after a certain 
time, and right before the resonant capture (at t « 2300) its in- 
ward migration is even slightly reversed, which is not seen in the 
figure. The action of the inner disk means th at we might expec t 
a slow outward migration of the inner planet (ICrida et al.l l2008). 
However, this does not play any role in the resonant capture in 
our case, as it is the outer planet that migrates towards the inner 
one. 
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Fig. 4. Dynamical behaviour of the two giant planets embed- 
ded in the surrounding protoplanetary disk: semi-major axes and 
eccentricities. During the first 500 orbital periods (of the inner 
planet), the planets are kept fixed to obtain a steady state in the 
disk. The upper panel shows the behaviour of the semi-major 
axes, while the bottom panel shows the time evolution of the 
eccentricities (red line corresponds to the inner, green line to 
the outer giant planet). The capture into the 3:1 MMR occurs 
around 2300 periods of the inner planet. The planetary masses 
(mi = 3.15My and m 2 = 1A6Mj) are taken from Fit II. 



If the effect of the inner disk is not taken into account, the in- 
ward migrating outer planet always cap tures the inner p lanet into 
an MMR. In this case, the results of lKlev etaf] d2004l) clearly 
demonstrate that the final system is either in 3:1 or in 2:1 MMR 
mainly depending on the speed of the outer planet's migration. 
A fast migration of the outer planet may result in its crossing the 
3:1 MMR without capture and ending the migration in the more 
robust 2:1 MMR. 

In the case of the large planetary masses, around 2300 pe- 
riods of the inner planet, the planets enter into the 3:1 MMR. 
The evolution of the semi-major axes and eccentricities during 
the migration and after the resonant capture is shown in Fig. [4] 
After the resonant capture, the eccentricities first increase and 
then tend to oscillate (around their mean values). However, the 
oscillations of the eccentricities remain bounded approximately 
in the range given by the numerical integration shown in Fig.Q] 
The inner disk certainly does contribute to the apparently very 
effective damping of the inner planet's eccentricity, as shown in 
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500 1000 1500 2000 2500 3000 3500 4000 4500 
time [period of the inner planet] 

Fig. 5. Dynamical behaviour of the two giant planets embedded 
in the surrounding protoplanetary disk for Fit II. The evolution 
of the resonant angle 6*3 and Aa> are shown in the top and bottom 
panels, respectively. After entering into the 3:1 MMR, the orbits 
of the giant planets exhibit apsidal corotation. 

the format ion scenarios of the reso nant systems HD 73526 and 
GJ 876 bv lSandoretaT] (120071) and lCrida et ail (l2008h . 

According to our simulation, the resonant angles (9,- (i = 
1 , 2, 3) and the corotation angle Aa> show libration, f?i librates 
around 0°, and the remaining angles around 180°, see Figure [3] 
The dynamical behaviour of our modelled system is qualitatively 
the same as the observed one, displayed in FigureQ] even though 
in our case the libration amplitudes are somewhat smaller. 

Interestingly though, the behaviou r of the resonant an gles 
does not match the theoretical results of Beauge et al. ( 2003]) and 
the numerical simulations done by iKlev et alJ (120041) . Although 
the eccentricity of the inner planet exceeds the limit of e\ ~ 0.13, 
the libration of #3 and Aa> is found to be around 180°. As al- 
ready mentioned, this might be the consequence of the system 
not yet reaching a stationary solution. In what follows, we study 
this unexpected result through numerical integration of the grav- 
itational three-body problem using properly parametrized non- 
conservative forces for migration. 

5. Toward the stationary solutions of the 3:1 MMR 
in HD 60532 

In the previous section we presented a reasonably well-working 
formation scenario for the system HD 60532. The eccentricities 



of the giant planets and resonant angles (including the corota- 
tion angle) of the formed system behave very similarly to the 
results of numerical integration using initial conditions from Fit 
II. In this section, we intend to explain why the obs erved and the 
forme d sy stems do not fit th e previous studies of Beauge et alJ 
d2003l) and lKlev et ail (j2004). To do so, we performed two series 
of simulations; in the first one, only the outer planet is forced to 
migrate inward, while the inner planet does not feel any addi- 
tional damping forces. In the second one, the inner planet feels 
the (accelerating) effect of the inner disk as well. 

The additional non-conservative force responsible for mi- 
gration and eccentrity damping can be parametrized by the mi- 
gration rate a/ a and the eccentricity damping rate e/e, or by 
the corresponding e-folding times t„ and t ? of the semimajor 
axis and eccentrici ty of the outer planet (see lLee & P eale 2002; 
iBeauge et alj |2006, for two possible approaches). I n our simula- 
tions, we used the force corrections as suggested bv lLee &' Peale 
(2002). The relations between the damping rates and e-folding 
times are a/ a = -l/r a , and similarly for the eccentricities 
e/e = -1/t c . One can define the ratio between the e-folding 
times K - t„/t,, (or e/e = -K\a/a\), which according to 
iLee & Peale! y002) determines the final state of the system in 
the case of a sufficiently slow migration. 

If the inner planet is influenced by an inner disk, it may be 
forced to migrate outward and its eccentricity is also damped. 
This effect can also be modelled by using a (repelling) non- 
conservative force parametrized by a/a (having a positive sign 
for outward migration), and e/e for the eccentricity damping. 
A detailed study of the resulting dynamical effec t of these 
param etrizations has been presented recently by ICrida et ail 
d2008l 

In the first series of simulations, only the outer planet mi- 
grated on different timescales t„, while its eccentricity was 
damped on timescales T e such that their corresponding ratios 
(K = r a /T e ) were always between K ~ 1 .. .10, which is typ - 
ically found in hydrodynamical simulations dKlev et al.l l2004). 
In the second series of runs, besides the migration of the outer 
planet, the inner planet's eccentricity was also damped, while its 
semi-major axis was not influenced. Typical cases of these sim- 
ulations can be seen in Figs.|6]and|7]ln these particular cases the 
following migration and eccentricity damping timescales were 
applied; for the outer planet 02/^2 = -5 X 10~ 5 year~', <»2/e2 = 
-2 x 10~ 4 year _1 (or T fl , = 2 x 10 4 , r e , = 5 x 10 3 years), and if the 
inner disk's influence on the inner planet was also taken into ac- 
count, we used for the inner planet h\/a\ = (as its semi-major 
axis was not influenced) and ei/ei = -2 x 10 _4 year _I . 

Comparing Figs. [6] and [7] it can be seen immediately that, 
in the first case when the inner planet was unaffected, the res- 
onant angles behave according to the theoretical predictions of 
Beauge et al. (2003) (asymmetric libration of the resonant an- 
gles), while in the second case the libration of the resonant an- 
gles is around 180° and 0°, respectively. The inclusion of the 
inner disk helps to keep the inner planet's eccentricity at lower 
values close to those coming from the observations. In Fig.|6]the 
eccentricity of the inner planet tends to 0.5, while in Fig.|7]it will 
show only small amplitude oscillations around the value 0.3. 

Finally, we can conclude that the effect of an inner disk in the 
three-body simulations (with dissipative forces) alters the final 
state of the resonant configuration (here a 3:1 MMR), keeping 
the resonant angles far from their mean libration, thus impeding 
the system for reach ing the equilibrium solutions suggested by 
Beauge et al. (2003). In both cases the systems will remain in 
their final configrations, when the damping effect of the disk is 
reduced and eventually non existent. 
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Fig. 6. Behaviour of the eccentricities (red for the inner, green 
for the outer planet), the resonant angle 6*3, and the corotation 
angle Aa> when only the outer planet migrates inward and the 
inner planet is not affected. 
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Fig. 7. Behaviour of the eccentricities (red for the inner, green for 
the outer planet), the resonant angle #3, and the corotation angle 
Aco when the outer planet experiences migration and eccentric- 
ity damping, and the inner planet only the eccentricity damping 
effect of an inner disk. 



6. Summary 

After the discovery of the system H P 60532 dDesort et alJ2008h . 
the thorough dynamical analysis of lLaskar & Corr eia ( 2009]) has 
reliably established the first resonant system that contains two 
giant planets in a higher order mean motion resonance, here 3:1. 
Until now, only one system, 55 Cancri, has been proposed as a 
candidate for a system where two giant planets might be in 3:1 
MMR. On the other hand, for 55 Cancri, the existence of the 



resonant solution could not be confirmed by iNaef et al] ( 20041). 
and a n ew non-resonant solution was later found by Fischer et al. 

Thus according to our present knowledge, HD 60532 is 
the only know n system containing gian t planets in the 3:1 MMR. 

Recently. iLask ar & Correia (2009) has improved the orbital 
solution found bv iDesort et al.l (12008). In their new fit that as- 
sumes a relatively small inclination between the common orbital 
plane of the planets and the tangent plane of the sky, i = 20°, 
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the giant planets have quite high masses: m\ = 3.15My and 
m 2 = 7.46M/. 

The main path to forming a resonant planetary system is 
thought to be by convergent migration of planets. For an alter- 
native scenari o, formation th r ough a scattering process has been 
proposed by iRavmond et alj (2008). Migration is the the result 
of dissipative forces originating in an ambient accretion disk, 
which act along with the mutual gravitational forces between the 
planets and the central star. Being a dissipative process, a suffi- 
ciently long-lasting migration process brings the system close to 
a stationary solution, which corresponds to a periodic solution 
of the system. The stationary (minimum energy) solutions in the 
systems yiel ding 3:1 MMR have b een studied by a Hamiltonian 
approach by Beau gi"et al.l (|2()03). It has been found that, de- 
pending on the ratio of the eccentricities, the resonant angles can 
exhibit antisymmetric (symmetric anti-aligned) or asymmetric 
libration. Thes e results have als o been confirmed by numerical 
simulations of lKlev etalJ J200l . 

In a first step, we integrated the system HD 60532 numeri- 
cally using the (observed) initial cond itions calculated from Fits 
I and II of lLaskar & Correial d2009h . We found that all reso- 
nant angles librate and that the system is in apsidal corotation. 
Since Aa> oscillates around 180°, observationally the system lies 
in an antisymmetric configuration. However, for the observed 
mean eccentricities, the Hamiltonian approach would suggest an 
asymmetric configuration. 

To model the formation of the system HD 60532 and find an 
explanation for this special antisymmetric configuration, we as- 
sumed that the system formed through a planet-disk interaction 
process. We performed a series of fully hydrodynamical simula- 
tions with low, intermediate, and high planetary masses /«,/ sin i, 
assuming i = 90°, 30°, and 20°, respectively. Since the eccen- 
tricity of the inner planet oscillates around a moderate mean 
value e\ ~ 0.3, we assumed that during the migration process 
an inner disk (between the inner planet and the central star) was 
present, providing an efficient damping mechanism on the in- 
ner planet's eccentricity. We found that the convergent migration 
through planet-disk interaction, which takes the planets into the 
3:1 MMR only occurred for the highest planetary masses. The 
dynamical behaviour of the resulting resonant planetary system 
is then indeed very similar to the one obtained from the radial 
velocity observations. Through our full hydrodynamical simula- 
tions, we supp ort the small inclination i = 20° of the system as 
suggested by Laska r & Correial d2009t) . 

To understand why the system does not appear to be in a 
miminum energy configuration to the 3:1 MMR, we performed 
a series of dedicated 3-body simulations with additonal forces 
emulating the effects of the protoplanetary disk. In particular, we 
studied situations with and without the influence of an inner disk. 
We found that it is exactly the effect of the inner disk that distigu- 
ishes between an antisymmetric or asymmetric final configura- 
tion. Here, its presence is responsible for the system not reaching 
the minimum energy asymmetric configuration. Nevertheless, 
we point out that the solutions obtained in the antisymmetric 
state are indeed stationary solutions, which are stable when the 
effects of the disk are reduced. 

In related earlier works, we had already established that the 
observed smallness of the mean eccentricities in resonant plane- 
tary systems can also be attributed to the effect of an inner disk 
(ISandor et al.l2007l:ICrida et al.l2008h . We conclude that an inner 
disk during the migration process has directly observable conse- 
quences on the post-formation dynamical behaviour of resonant 
planetary systems in mean motion resonances. 
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